The Jacobi-Rosochatius problem on an ellipsoid: the Lax 
representations and billiards 



BOZIDAR JOVANOVIC 



Abstract. The Lax representations of the geodesic flow, the Jacobi-Roso- 
chatius problem and its perturbations by means of separable polynomial po- 
tentials, on an ellipsoid are constructed. We prove complete integrability in 
the case of a generic symmetric ellipsoid and describe analogous systems on 
complex projective spaces. Also, we consider billiards within an ellipsoid under 
the influence of the Hook and Rosochatius potentials between the impacts. A 
geometric interpretation of the integrability analogous to the classical Chasles 
and Poncelet theorems is given. 
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1. Introduction 

A well known Jacobi's problem describes the motion of a material point on a 
n-dimensional ellipsoid i?" 

(1.1) = {xe R"+i I {A-^x, x) = 1} 

under the influence of the Hook elastic force —ax. Here A = diag(ao, ai, . . . , a„) is 
a positive definite matrix and cr is a real parameter [24L 136) . The Lagrangian of the 
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system is L{x, x) = ^{x, x) — ^{x, x) . The corresponding Euler-Lagrange equations 
are 

(1.2) x = \A^^x-ax, 

where the Lagrange multipher is A = —{{A^^x^x) — a) / {A^"^ x , x) . 

By introducing the momenta y = dL/dx — x, we can write (|1.2p in the Hamil- 
tonian form 

(1,3, ,^y. i.-™.^.-..-,. 

on the cotangent bundle T*E^ reaUzed by the constraints 

(1.4) Fi = {A-^x,x) -l = Q, F2 = {A-^x,y) ^0. 

The equations p.3p are Hamihonian with respect to the Hamihonian 

(1-5) H{x,y) = {y,x) - L(x,x)\.i=y = ^{y,y) + ^{x,x) 

and the Dirac-Poisson bracket defined by 

{F,,h}{F2j2}-{F2,h}{F,j2} 



(1.6) {/i,/2}d-{/i,/2}- 



{Fi,F2} 



where {•, •} is the canonical Poisson bracket on ]R^"+^(x, y) (see [3l 136) ). The 
Dirac-Poisson bracket of the coordinate functions are: 

•^iHj '^jVi r T e XiXj 



{x,,Xj}d = 0, {y^,yJ}D = rjz:^ r, {x^,yJ}D = 7-77:^ r- 

aiaj{A ^x,x) aiaj{A ^x,x) 

Note that the Hamiltonian flow and the Dirac-Poisson structure is defined not 
only on T*E'^ but on the whole M2n+2 without a; = 0. 

If we set (7 = 0, the system represents the geodesic flow on the ellipsoid (|l.ip . 
The geodesic flow and the Jacobi problem are one of the basic classical models 
of integrable systems [241 [3]. For the ellipsoids with distinct semi-axes, they are 
separable in the Jacobi elliptic coordinates (Aq = 0, Ai, . . . , A„) defined as follows 

[S [Ml nil [50]. Through every point x e R"+\ xq ■ xi a;„ 7^ 0, pass exactly 

n -f 1 mutually orthogonal confocal quadrics Qaq , Q\i , • ■ • , Qa„ (Aq < aq < Ai < 
ai < • • • < A„ < a„) given by the equation 



(1.7) Qa: ((A-A)-ix,x)=^^ = l. 

The original coordinates are defined up to a sign: 

In his celebrated paper |36j . Moser found the Lax representation L = [L,B\^ 
L = (^I„+i - ^) {A-x<E>x) (^I„+i - ^) , B = A-'x A A'^y, 
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of the associated system 

(1.8) x^^, y = <f = ~{y,A-^y){l-{x,A-^x))~{x,A-^y)\ 

The set $ = describes the set of tangents of the elhpsoid E" — Qq. Restricted 
to $ = 0, the system (jl.Sp has the following important property: the moving line 
p{t) = {x{t) + sy{t) I s € R} has a point of contact ^ with i?" along nonparameter- 
ized geodesic ^(r). Moreover, the eigenvalues rji, . . . ,?7„ of L (different from zero), 
define n quadrics , . . . , from the confocal family p.7p . such that p{t) is si- 
multaneously tangent to Qjj^ , . . . , (a variant of the classical Chasles's theorem 

El)- 

Although Moser's method is applied in the construction of various integrable 
models admitting Lax representations [1], Lax representations with a spectral pa- 
rameter of the geodesic flow and the Jacobi problem on £■", according to author's 
knowledge, are not given yet. 

In Section 2 we define double Jacobi flows and complex double Jacobi flows 
and construct two different Lax representations, a "small" 2x2 one fTheorem l2.ip 
and a "big" (n + l) x (n-|- 1) one fTheorem l2.2p . An appropriate restrictions of the 
double Jacobi flows and the complex double Jacobi flows are the Jacobi problems 
on the ellipsoids p.ip and 

(1.9) £;2n+l ^ 1^ ^ I ^^-1^^ ^ 

respectively, leading to the Lax representations of the Jacobi problem (Theorem 
13.11 Section 3). Note that a small 2x2 Lax representation has the usual Lax matrix 
form of the Jacobi-Mamford system (see |33L 145)1. 

The Jacobi problem on the ellipsoid (|1.9p is invariant with respect to the stan- 
dard T""'"^-action on C""''^Q It is well known (e.g., see [36, 28j), that the reduced 
flow can be naturally considered as a system describing the motion of a material 
point on the ellipsoid (jl.ip under the influence of the Hook and the Rosochatius 
potentials |40j 

1 " 2 

(1.10) V{x) = ^{x,x) + -J2^ 

(Section 4). The Lax representation is also invariant with respect to the standard 
T"+^-action and induce a small 2 x 2-Lax representation with a spectral parameter 
for the Jacobi- Rosochatius system on the ellipsoid E" (Theorem 14. ip . Note that 
if instead of the T"+^-reduction, we perform a S'^-reduction, with diagonally 
embedded in T""^^, we obtain a natural mechanical system on a complex projec- 
tive space (Proposition 14. ip . providing a class of examples of Hermite-Liouville 
manifolds (see Remark |4T|) [2gll44] . 

If all semi-axes of the ellipsoid (jl.ip are distinct, the Jacobi- Rosochatius sys- 
tem is separable in elliptic coordinates and has n independent commuting integrals 



^We use the complex notation to simplify the reduction procedure. Equivalently one 
can consider the real space R-^"+^ and the ellipsoid iJ^n+l fjggng(j jjy ^j^g matrix A = 
(ao, ao,ai,ai, . . . , an,an). 
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which are quadratic in raomenta. The geodesic flows on symmetric ellipsoids are 
studied in details in |10| . In general, the geodesic flows are integrable in a non- 
commutative sense. In Section 5, we prove complete integrability of the symmetric 
Jacobi-Rosochatius system both in a noncommutative (Theorem 15. 2p and the Li- 
ouville sense (Theorem 15. 3p . 

In Section 6 we consider the billiard system within an ellipsoid under the in- 
fluence of the potential (jl.lOp between the impacts. By a slight modification of the 
results given by Fedorov |21j , we describe the billiard mapping and the Lax repre- 
sentation fTheorem lG.ip . A geometric interpretation of the integrability analogous 
to classical Chasles and Poncelet theorems is given (Theorem 16. 2p . 

Finally, in the last section, by using the 2x2 Lax representation, we give the 
hierarchy of the Lax representations for the separable potential perturbations of 
the Jacobi-Rosochatius system on the ellipsoid E". 

2. Double Jacobi flows 

Together with the Jacobi problem, let us consider a system defined by the 
Lagrangian L(x,^,i:,^) = (i,^) — a{x,^) and constrained on the hypersurface 

S = {{X,0 E M2n+2 I (^^^-1^) -1=0}. 

The Euler-Lagrange equations are 

X = XA^^X — (TX, 

where the Lagrange multiplier is 

As in the case of the Jacobi problem, we can write the corresponding Hamilton- 
ian equations on T*E. Canonically conjugate momenta to {x, ^) are r] = dL/dx = ^, 
y ~ dL/d^ ~ X. The Hamiltonian function of the system is given by the Legendre 
transformation: 

H(x, ^, T], y) = {x, 77) + {i,y) - L = {y, nf) + (j{x, £,), 
and the Hamiltonian equations are 

(2.1) x = y, y = -— ^ — - — A x - ax, 



{A-^x,0 

(2.2) ( = ^ = _____A 



The symplectic structure on 

(2.3) T*E : Gi = {x, A~^0 -1 = 0, G2 = (y, A'^O + i^, ^-^v) = 

is the restriction of the standard symplectic structure on R^"+''(a;, ^, 77, y). The 
corresponding Dirac-Poisson structure is given by (jl.6l) , where we should replace the 
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constraint functions Fi,F2 by Gi,G2 and {•,•} is the canonical Poisson structure 
on R^''+^{x,^,r],y). 

By the analogy with the double Neumann system (see |41j ). we refer to (j2.1l) . 
p.2p as a double Jacobi flow. 

We can also consider the complexified phase space, an affine 2n + 2-dimensional 
variety in C^^+^ix, C, V, v) defined by (gS]). Then we refer to ([23]), as a com- 
plex double Jacobi flow (here, the time parameter t can be both: real or complex). 

Let A\ — diag(A — aq, A — ai, . . . , A — a„) and 

n ^ 

q,ix,0^{A-,'x,O^T.Y^- 

2 = 

Theorem 2.1. The (complex) double Jacobi flow (j2.ip . (j2.2D implies the matrix 
equation 

(2.4) C{X) = [C{X),A{X)] 
with 2x2 matrices depending on the parameter A 

-1 - q\{x,£,) -q\(y,0 

Proof. To simplify the calculations, consider the time reparametrization 

(2.5) dt = {A~^x,^)dT. 

The double Jacobi flow in the new time gets the symmetric form 

x' = {A-''x,Oy, 

y' ^{(7- {A-^y, T^))A-^x - a{A-^x, Ox, 

(2.6) e = {A-^x,ov, 

V'^ia- {A-'y, 7^))A-'^ - a{A-'x, 00 

where (.)' ^ ^(.) ^ {A-^x,Oii-)- 

By using the constraints (j2.3p and the identity 

A-'A^' = A^'A-' = (A-' + A^')/X, 

we obtain the relations 

{qx{x,v)y = {q,{y,Oy - ^(1 + {A-'x,0){<y-{A-'y,v)) 

+ {A-^x,0{{y,A^'r,) -a{x,A^'Ol 

{qx{y,v)y = \{{A-,'x,v) + {A-,'y,0){<y {A-'y,v)) 

-a{A-^x,0{{x,A-'7^) + {A-'y,0), 
{qx{x,Oy = (A-2x,e)((x, + {A-^'y,0), 

Now, the reader can simply verify C — [C, {A^'^x,OA]. □ 
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Therefore, the coefficients of the invariant polynomials of the matrix a(X)C{X) 
are the first integrals of the system. If all of Ui are distinct, the integrals can be 
written in the form 

gi ^ VtCi - XiT]i, i 0,l,...,n. 

Theorem 2.2. The (complex) double Jacobi flow (12. ip . (|2.2p restricted to the 
invariant variety 

(2.7) {A-'x,r])^{A-'y,O=^0. 
implies the matrix equation 

(2.8) /:*{X) = [A*iX),C*W] 

with {n + 1) X (n + 1) matrices depending on the parameter X 

£* {X) ^ X{y ^ - X 1]) + y (S) Tj + ax ^ - (tA - A^A, 

^*(A) = I (A-^y ^ A-^C - A-^x <S) A'^i] + XA'^) . 
{A-'^x,Q 

Proof. From (|2.6p we get 

ix(^Cy = {A-'^x,0{y<»C + x(^T]), 

(y^vY = i<^- {A-^y,r])){A-^x rj + y <g) A-^() - a{A-'^x,Oix <E) rj + y C), 
iy(^^- xi^rj)' = (a - {A^^y, 7j)){A-^x (g) ^ - x A^^^). 



Thus, 

(£*(A))' = (a - {A-^y, 77)) {X{A-'x ®^^x® A-^0 + ^'^^ ^V + V® ^"'C) • 

Further, let us denote 

n = A^^y (g) A^^S, - A^'^x (g) A^^ri. 

Then, by using the constraints (12.31) . we get 

[ft, y®rf\ = {A^^y, ClA~'^y <»V + {A~'^x, ri)y ® A^^r) 

-{A-^-q, y){A-^x ®-q + y® A^^i), 
[rj, X (g = A^^y (^"S, x)A^^x {A^^y, Cix g) A^^S^ +x® A^^rj, 

[ri, y g) ^ - X g) 77] ^ y (g> A^'^f] - A^'^y (g) T] + {A^^x, rj) {A~^x (g 77 - a; (g A^^yy) 
+ {A-^y, {A-^y (g ^ - y (g A^^^ - (A^^y, ry) (A-^x (g ^ - x g) A-^^), 
[ft, A] ^ A^^y g) C - A'^^x (E)ri-y(E) A^'^^ + x (g) A^'^rj, 
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implying 

[{A-'^x, e)A*{\), C*] = [rL + \A-^,X{y ®i,-x®ri)+y®j^ + ax®S,-aA- X^A] 

= (o- - {A^^y, ri)) {\{A^^x ®£^~x® A^^£) + A^^x ® rj + y ® A^'^C) 

+ {A^'^y,Cl{A^^y®'n ~ ax® A^'^S) + {A^^x,rj){y ® A^^-q - aA^^x®^) 

+\ {{A-'^x,ri){A-'^x ® -q - X ® A-'^f]) + {A-^y,C){A-'^y ® ^ - y ® A'^Cl) . 

It remains to note that (j2.7l) defines an invariant manifold of the system. In- 
deed, we have 

{{A-'x, ry))' - (A-^x, {y, A-'r^) + [a - {A-'y, ^MA'^^, - '^(^"'^, = 0. 

□ 

Note that the matrix C{X) is invariant under the transformations of the phase 
space given by the (R*)"+^-action (respectively (C*)"+^-action): 

while the matrix £*(A) is K*-invariant (respectively C*-invariant): 

[xiXiiVi^Vi) ^ — > {sxi,s^^£,i,s~'^rii,syi), s ^ 0, i = 0, 1, . . . , 

3. The Lax representations of the Jacobi problem 

The equations 

x = £., y = v 

define the invariant manifold of ()2.ip . (j2.2p . so the double Jacobi flow contains as 
a subsystem the Jacobi problem (jl.3l) . Also, for a; = ^, 7/ = 77, the condition (|2.7p is 
satisfied and the above theorems imply Lax representations for the Jacobi problem. 
In particular, when we set cr = 0, we get the Lax representations for the geodesic 
flow on a ellipsoid. 

Theorem 3.1. The Jacobi problem (jl.3D implies the matrix equations 

(3.1) C{X) = [C{X),A{X)] 
and 

(3.2) C*iX) ^[A*{X),C*iX)] 

with 2x2 and (77 + 1) x (n + 1) matrices depending on the parameter X 
q\{x,y) qx{y,y)+a^ 



CiX) = 



-l-qxix,x) -qx{y,x) 



{{a-{A-^y,y))~a{A-^x,x} 
■^^'^{A-^x,x)[{A-^x,x) 

£*(A) ^ X{y A x) + y (g) y + ax (g) X - (jA - X'^A, 



respectively. 
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From the Lax representations we obtain the well known form of the integrals 
of the Jacobi problem given by Moser |36; . Note that the 2x2 Lax matrix £(A) 
has the usual Lax matrix form of the Jacobi-Mamford systems (see |33L 145) 1. Also 
note, if all are distinct, the Lax representations (13.11) and (13. 2p are equivalent 
to the Jacobi problem (|1.3p . up to the action of groups ^2^^ , Z2 generated by 
reflections 

{xi,yi)i — > {siXi,Siyi), = ±1, z = 0,l,...,n, 
{x,y)i — > s{x,y), s==±l. 

The time reperemetrezation (|2.5p . for x — S,, coincides with the time reparametriza- 
tion dt = Xi . . . Xndr used in the integration of the geodesic flow [48,. ^ 20, . 



4. The Jacobi-Rosochatius problem on an ellipsoid 



Consider the complex double Jacobi flow (12. IL (12.21) on the invariant real va- 
riety 

(4.1) x = z, £, = z, y=p, r] = p. 

The equations 1^^, dHJ) define the cotangent bundle r*£;2n+i ^^le 2n + 1- 
dimensional real ellipsoid in the complex space (jl.9l) 

(4.2) T*^2n+i. (A-iz,z)-l = 0, {A-^z,p) + {A-H,p) ^Q. 



+ 77'2n— 1 



The complex double Jacobi flow restricted to T*E 

(4.3) p= -^^^^^^^-^A-^z-az, 

describes the motion of a material point on i?2"+^, under the influence of the 
elastic force —az. From Theorem 12.11 we see that the system (|4.3p implies the 
matrix equation 

(4.4) £(A) = [/:(A),^(A)], 
where 



C{\) = 



q\{z,p) qx{p,p)+a 
-I - q\{z,z) ~qx{p,z) 

1 / j-{<7-{A-^p,p))-a{A-^z,z) 



A ^ 

-2. 



{A-^z,z) \{A^^z,z) 



The system (|4.3p . as well as the Lax representation (14. 4p . is invariant with 
respect to the Hamiltonian torus action on x*E'^"^^: 

(4.5) {zk,Pk)' — > e"^''{zk,Pk), fc = 0, l,...,n, 
with the momentum mapping 

z z 

(4.6) h.= {ho,hi,...,hn), hk = --gk = -^{zkPk - PkZk)- 



THE JACOBI-ROSOCHATIUS PROBLEM ON AN ELLIPSOID 



9 



4.1. The Jacobi system on a complex projective space. In particular, 
the Hamiltonian h = Hq + ■ ■ ■ + hn induces a S'^-action 

(4.7) iz,p) ^ e'^{z,p). 

The syniplectic reduced space h~^{Q)/S^ is simplectomorphic to the cotangent 
bundle of the complex projective space P". The reduced system is a natural me- 
chanical system with the kinetic energy determined by the ellipsoidal metric on 
P", the submersion metric with respect to S^-action, under the influence of the 
"elastic" force. From Theorems [O and [El where we identify C"+i M2"+2 
set A = diag(ao, oq, ai, ai, . . . , a„, a„), we have integrability of the Jacobi problem 
(|4.2p for an arbitrary choice of Oj. Similarly, the reduced flow is integrable. A 
general treatment of the integrability of the reduced systems is given in |52L I25j . 
Here, the orbit of the S'^-action are tangent to the isotropic tori 7", which lay in 
/i~^(0), and T/S^ are invariant isotropic tori for the reduced flow. 

In the next statement we describe the reduced system on P". Let tt : C""''^ \ 
{0} — >■ P" be the canonical projection oi w = {wq, . . . , Wn) G C""*"^ \ {0} to [w] = 
[wq Wn] G P", with respect to the C*-action. 

Proposition 4.1. The reduced Jacobi problem on P" is a natural mechanical 
system with the kinetic energy determined by the metric 

n X. XW {w,Aw){X,AX)-{X,Aw){w,AX) 

gA(TT*X,TT^A)U^] ^ — — , A e J^(L \ {Oj 

[w, Aw){w,w) 

under the influence of the potential field V{[w]) — <y{w, Aw) /2(w, w) . 

Proof. Under the change of variables Zi — yfalwi, the Jacobi problem trans- 
forms to the system on a sphere 5'2"+i : {yj^ y}) = i with the Lagrangian function 

L{w,w) ~ -{w,Aw) — —{w,Aw). 

It is well known that the reduced system on h^^{0) / = T*P" can be described 
by the S'^-La grange- Routh reduction with a zero value of the S^-momentum map- 
ping (e.g., see [3]). Since it is convenient to work with homogeneous coordinates, 
we firstly extend L to the C*-invariant Lagrangian function 

(WjAw) a{w,Aw) 



(4.8) L 



2{w,w) 2{w,w) 



defined on C"+^ \{0} and then perform the Lagrange- Routh reduction with respect 
to the C*-action. 

Obviously, the reduced potential is T^(['i«]) = (t{w, Aw) /2{w,w). Further, the 
Lagrangian (|4.8p defines the Riemannian and Hermitian metrics 

2{w,w) {w,w) 

respectively. The reduced system is a natural mechanical system {F"^ , gA,V{[w])), 
where gA is the submersion metric with respect to the C*-action [3]. 
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At every w G C"+^ \ {0} we have a decomposition 

T„C"+i \ {0} = C"+i = ® 

where Vw is the tangent space of the orbit of C*-action through w {vertical space 
at w) and is its ij^-orthogonal complement {horizontal space at w). Since Vw is 
the complex line through w, its real g^^-orthogonal complement coincides with its 
/lA-Hermitian orthogonal complement 

n^ = {x\hA{x,w)^Q}. 

Let X,Y ^ T„,C"+^ \{0}. By definition, the submersion metric gA(7r*X, 7r*)|[^] 
on P" by is equal to gA{X' ,Y')\w, where X' and Y' are horizontal components of 
X and Y 



hA_{w,w) ' hA{w,uu) 



Therefore 



gA(7r,X,7r,X)|[,,] = gA{X\X')\^ = hA{X\X')\.^ - hA{X\X) 

u (Y Y\ hA{w,X)hA{X, w) 

= hA[X,X) — 

hA{w,w) 

{w, Aid) {X, AX) - {X, Aw) {w, AX) 
{w, Aw){w^ w) 

□ 

Remark 4.1. Note that (P", qa) is an example of a Hermite-Liouville manifold 
(see |23L I44j). In particular, if A is the identity matrix, then gA is the standard 
Fubini-Study metric on P". The integrability of the geodesic flow of the Fubini- 
Study metric is proved by Thimm |43| and Boyer, Kalnins and Winternitz [7]. 
Further, besides the Hook potential in (|4.3|) . we can add other separable polynomial 
potentials V{z) (see the last section), and by the Maupertuis principle (e.g., see [3]), 
(P", (c — V{[w]))gA), c > maxy([w]) are examples of Hermite-Liouville manifolds 
as well. 

4.2. Reduction to the Jacobi-Rosochatius problem. Now we shall per- 
form the reduction with respect to the torus action (|4.5p . It is well known that we 
obtain a natural mechanical system under the influence of the Rosochatius potential 
[40ll36l[28] . 

Introduce the canonical polar coordinates {xk,^k,yk,hk), where hk are given 
by and (pk,yk by 

(4.9) Zfc = Xfee''^^ pk = y^e''^'^ + i^^e'^" , fc = 0,l,...,n. 

Xk 

The Hamiltonian of the Jacobi problem in new coordinates reads 



k 

,2 ■ 

-n -^k 



k=Q 
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(4.10) Xk = yk: Vk = T-r^ ^ flfe Xk - (TXk 



We see that ipk are cyclic variables of the system. Consider the level set of the 
momentum mapping 

(T*i?2»+i)^, : h=(A.o,Mi,---,Mn) 

(some of constants /i^ may be equal to zero). 

The equations of the system on (r*^'^""''^)^ separate on reconstruction equa- 
tions 

Xfc(y9fc=^fc, fc = 0, l,...,n. 
and the reduced system on (T*£'2"+^)p/T"+^ in variables {x,y): 

{A-'^x,x) 

where fi/x = (hq/xq, ^ii/xi, . . . ,//„/x„). 

Note that, according to the definition (|4.9p . we have Xk > 0. Also, the variables 
{x,y) satisfy the constraints (|1.4p . However, we can consider the system (|4.10l) on 
the whole cotangent bundle (jl.4p . Then it represents the motion of a material point 
on a ellipsoid with the potential energy V{x) having two terms: the Hook and the 
Rosochatius potential (|1.10p . 

Applying the change of variables (|4.9p to the Lax representation (|4.4p on the 
invariant set {T*E'^"^-^)^, after subtracting the multiple of the identity matrix from 
£(A), we get the following statement. 

Theorem 4.1. Suppose that the eigenvalues at of the matrix A are distinct. 
Up to the action of the group 11^^^ generated by the reflections 

[xi,yi)^ — > {siXi.Siyi), s.^ = ±1, i = 0, 

the J acohi- Rosochatius problem (|4.10p is equivalent to the matrix equation 

(4.11) t{\)^[C{\),A{\)] 

with 2x2 matrices depending on the parameter A 

C{X)=( l^(^^y) 9A(z/,y) + gA(^,^)+CTA 
- q\{x,x) ~qx{y,x) y' 

i f i(a-(A-iy,y)-(A-if,f))-a(A-2x,x) 



■^^^^ ^ {A-^x,x) {{A-^x,x) 
The integrals obtained by (j4.1ip can be written in the form 

(4.12) /, = yf + + 1 + E ^ (iy^-^ - y.--^f + 4? + ^ V 

i = 0, 1, . . . , n. This is a commutative set of functions, both with respect to the 
canonical Poisson brackets {•, •} and the Dirac-Poisson bracket {•, - Id- 

Let us note that a possible alternative approach in the construction of the 
Lax representations for the Jacobi-Rosochatius problem (|4.10p is by using the Lax 
representations of the Neumann problem (e.g., see j41| ) and the well known cor- 
respondence between the Neumann problem and the geodesic flow on an ellipsoid 
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[311 128) . An algebro-geometric study of the Neumann system on a sphere S" with 
the addition of the Rosochatius potential is given in [18 . 

5. Symmetric ellipsoids. Complete integrability 

Consider the case of a symmetric ellipsoid: 
aj=ao, i e /o = {0, . . . ,fco - 1}, 
(5.1) flj^ai, i e /i = {fco, . . . , fco + fci - 1}, 

tti = ar, i <E Ir — {ko + • • • + fcr-l, ■ ■ ■ , ko + ■ ■ ■ + kr-l + kr ~ 1}, 

ai aj,i j, k() + ■ ■ ■ + kr = n + 1. Then the Jacobi- Rosochatius problem (j4.10p 
is not equivalent to the Lax representation (|4.1ip . but still implies it. From the 
Lax representation (|4.11l) . the invariants of £(A) are integrals of the flow. By using 
the relation 

detCiX) = {l + qxix,x))(qxiy,y)+qx(-,-) +C7) -qx{x,yy 
we get the integrals 

(5.4) Ps= E P-i^ k,^\h\>2, 

i,jels,i<j 

where are given by 

22 22 

Pij — [ViXj — XiUj) H "2 I -2—, 

while for \Is\ = 1 we set = 0. Whence, we have p nontrivial integrals among Pg, 
where p is the number of sets Is for which ks = \Is\ > 2. 

In terms of integrals fg , the Hamiltonian of the system can be express as 



1 ~ 



H 

2 

s=0 

Theorem 5.1. Apart of the integrals arising from the Lax representation, the 
rational functions 

(5.5) -Ps,ij '■— Pij J *jj S Is 

are integrals of the Jacobi- Rosochatius problem (j4.10p . The functions /s,Ps are 
central functions within the set of integrals J- = {/s, PsSj}- 

{4,/.Jl3 = 0, {4,P,Ji3=0, {Ps„Ps,}d^G, 

{/.i,^.2,u}l3=0, {Ps,,Ps,,.j)d=Q. 
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Also, the functions Psi.iiji o,nd Ps2.i2j2 mutually commute for distinct si and S2, 

(5.6) {PsiMJi^Ps2a2j2}D ^0- 

Proof. The theorem can be verified by a straightforward calculations. In- 
stead, we consider the Jacobi-Rosochatius problems on an one-parametric family 
of deformed, non-symmetric ellipsoids 

El': {A-^x,x)^l, A, =diag(4,...,a^), < 7^ a^', i ^ j, 

where 

lim al = Oi, 

e-i-O 

and al are smooth functions defined on some interval (—A, A). 

Let be the associated Dirac-Poisson bracket (|1.6I) and let // be the 

corresponding integrals (|4.12l) . 

(5.7) {/^,/;}^ = 0. 
Define 

/:-E/^-Efj/?+-? + ^ + E^V . = 0,...,., 

Pl^o={<-a))ft, 

From ([STT]) we get {fl^Jl^Yo = £• Since lime^o ft = /s, by taking the 

limit we obtain 

{fsi,fs2}D = 0. 

On the other hand, limj^o Pt ij '^^'^ be singular and depends on the deforma- 
tion. Suppose that 



a- — a, 

lim^ 7 = 0, l^i,j- 

e->o a\ - a\ 

Then lime^o = PsSj and from {fl^,Pl^.i.jYD = fo^' all e, we get 

{/si,^S2,y}-D = 0. 

Therefore 

{fs,.Ps2}D= E {/.i,^..,.j}l?=0 
i.,jels2 ,i<j 

and 

s=0 

that is Ps,ij are integrals of the system. Similarly, for si 7^ S2, G /s^, 12, J2 € 
Is2 , we can always find a perturbation such that 

lim ~ °f ^ 0, Zd^i<i,jd, (i=l,2. 
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Therefore, linif^o Pg^Aajd ^ ^Sd.^djd^ d= 1,2, and we have 

{-fsi 1 ^S2,i2 J2 } ~ ^ ] Ps2,i2j2}D — 

il,jielsi,il<jl 

{Psi,Ps2}^ X! {Psi,Ps2,12]2}d ^ 0- 

It remains to prove {Ps, Ps,ij}D = 0, |/s| > 3. For shiiphcity, assume s = 0, i = 
0,j = 1. Consider a deformation having the property 



(5.8) 



hm°" 


-«! 


= 0, 


e-s-o af) 










Km — 




= 0, 


e->0 a| 




Km — 




= 1. 




-«! 





For example, we can take Oq = aoi = o:q + f-^, — ao + + e^, a| = ao + e"^ + 
+ e, . . . , ttko = ao + + €^ + (fco — 2)e. Subsequently, we get 

lim^o.o.i = lim(ao - a'i)/o' = Po,o,i> 

e— >-0 e— fO 

1™ ^0 2 = lim,(a2 - ao)/2 = Po,2,0 + ^0,2,1 = -Po,0,2 + Po,l,2 

and, consequently, 

(5.9) {Pofi,l,Po,l,2 + Po,0,2}D = 0. 

Next, we take = ao, af = ctQ + e^, al = Q^o + + a"^ = cto + 2e^ + e, 
a| = ao + 26^ + 2e . . . , a^, = ao + 26^ + (/cq - 2)e. Then 

l"^,4^-0, /^0,1, lim^4^^0, 1^2,3, 

and lime^o ^"0,0,1 = -Po^o.i, lim<:_>o ■^0,2, 3 = ^0,2,3- Thus, we get 

(5.10) {Pofi,i,Po,23}D^O. 

Finally, repeating the arguments given for (|5.9p and (|5.10l) , we obtain the com- 
muting relations 

(5.11) {Ps^^J,Ps^^k+Ps^Jk}D = 0, {Ps,,j , Psm} D = , hj^k,l, 

implying 

{Ps,^j,Ps}D = {Ps,^3,Ps,kl}D = J2 {P'>,^3 ' D + {Ps,^3 , Ps,jk} D = 0. 

k,l£ls,k<l k£ls 

□ 
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In particular, for the case of the Jacobi problem (where we set fJ^o = ■ ■ ■ = i-J-n = 
0), or for the case of a geodesic flow on a symmetric ellipsoid (/io = • • • = /i„ = cr = 
0), the system is invariant with respect to the 5'0(fco) x • • • x SO{kr) action and 
the integrals (|5.5|) reduce to the squares of the Noether integrals 

(5.12) ^s,i3 ^ ViXj ~ Xiyj, i < j, ijels, 
while the central functions Pg reduce to the invariants 

(5.13) "^hr 

i,jels,i<j 

The detail study of the geodesic flow is given in [10' (the case n = 3 can be 
found in 1111 ) . It is proved that among the central functions fs, Ps, with /xq = 
• • • = Mri = = 0, there are r + po independent ones, while that among fg and the 
Noether integrals (I5.12[) there are 2n — r — p independent ones. Therefore, for a 
sufficiently small parameters pi and a, we get that the dimensions of linear spaces 

F{x,y) = {Xj:^{x,y),Xp^ .^{x,y) \s = 0,1, ... ,r, i,j e Is), 
K{x, y) = (Xj^ (x, y), Xp^ (x, y) | s = 0, 1, . . . , r) 

are at least 2n — r~p and r + p, respectively, at a generic point (x, y) G T* E^^. Here 
Xf denotes the Hamiltonian vector field with respect to the Dirac-Poisson bracket 
(|1.6p . According to (|5.2I) we have the relation 

s=0 s=0 ^ i=0 ' 

so dim K{x, y) < r + p. 

Since all object are rational functions, dim F{x, y) > 2n—r—p and dim K(x, y) = 
r + p, for a generic values of pi, a, {x,y) G T*E". As a result, we conclude that 
T = {fs, Ps,ij} is a complete set of integrals (a generic dimension of F{x, y) equals 
2n ~ r ^ p and dmiF{x,y) + dimK(x,y) — dimT*i?") and we can apply the 
Nekhoroshev-Mishchenko-Fomenko theorem on noncommutative integrability (see 
[3811311 [3]). 

Theorem 5.2. The Jacobi-Rosochatius problem (j4.10p (the Jacobi problem 
(II. 3p ) on a symmetric ellipsoid (|l.ip . (|5.ip is completely integrable in a non- 
commutative sense by means of integrals (|5.3p and (j5.5p (where we set pi — 0). 
Generic trajectories take place over r + p-dimensional invariant isotropic tori, 
spanned by the Hamiltonian vector fields Xj ,Xp_^. 

In general, noncommutative integrability implies the usual Liouville integrabil- 
ity by means of smooth commuting integrals [6] . Here, the integrals can be chosen 
to be linear functions of non-commuting integrals. 

Theorem 5.3. The Jacobi-Rosochatius problem (|4.10p on a symmetric ellipsoid 
(|l.ip , (|5.ip is Liouville integrable by means of integrals (15. 3p and 

(5.15) = Y -^"'^J' k = l,...,ks-l, s = 0,...,r, 
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where 

10, k = {0,...,fc}, 

11, k = {ko,...,ko + k}, 

Ir,k ^ {ko-\ 1- fcr-1, • . ■ , fco H 1- fcr-1 + HI- 
PROOF. According to (|5.6p . we have {Lg^^ki, Ls2,k2}D — for si 7^ S2, while 
from (|5.1ip we get 

{Ps,ij,Ls,k}D = 0, i,j e /s,fe, 

and, in particular, 

{is,fei, is.fcalr) = 

(note that Pg — Ls^k^-i)- 

It is clear that the integrals (j5.15l) arc mutually independent and their total 
number is 

{ka-l) + --- + {kr-l)^{n + l)-{r + l)^n-r. 
On the other hand, from the completeness of = {fs,Ps.ij}, f functions among 
fs are independent from the integrals Ps.ij (we have the relation (|5.14p among 
them). Therefore, the set of integrals {fs,Ls.k} is a complete commutative set on 
T*E'^. □ 

A choice of commuting integrals is not unique and (j5.15p is motivated by a chain 
of subalgebras method in the construction of commutative functions on Lie algebras 
(e.g., see Thimm [43]). Another complete families of commuting integrals can be 
obtained, for example, by using separable variables related to the degeneration of 
the elliptic coordinates (see [71 126p. 

6. Billiards inside ellipsoids 

6.1. Billiards: continuous and discrete description. Let {Q,g) be a n- 
dimensional Riemannian manifold and let D C Q be a domain with a (smooth) 
boundary F. Let tt : T*Q Q he a natural projection and let g^^ be the 
contravariant metric on the cotangent bundle. Consider the reflection mapping 
r : TT^^T — > TT^^F, y_ y+, which associates the covector y+ e T*Q, x e F to 
a covector y_ E T*Q such that the following conditions hold: 

|y+l==|y-l, y+-y--Lr. 

A billiard in D is a dynamical system with the phase space T*D whose tra- 
jectories are geodesies given by the Hamiltonian equations with the Hamiltonian 
iJ(x,y) = ^5~"'^(y,y), reflected at points x G F according to the billiard law: 
r(y_) = y+. Here y_ and y+ denote the momenta before and after the reflec- 
tion. If some potential force field V{x.) is added than the system is described 
with the same reflection law and Hamiltonian equations with the Hamiltonian 
if(x,y) = i5-i(y,y) + y(x). 

A function f : T* Q M. is a.n integral of the billiard system if it commutes with 
the Hamiltonian {{f,H} = 0) and does not change under the reflection (/(x, y) = 
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/(x, r(y)), X e r). The billiard is completely integrable in the sense of Birkhoff it 
it has n independent integrals polynomial in the momenta, which are in involution 
(see [29] 'I. 

It is well known that the billiard system within an ellipsoid E"~^ C M" under 
the influence of an arbitrary potential separable in elliptic coordinates is completely 
integrable in the sense of Birkhoff |29l I13LI14) . Moreover, the manifold T*D/r \ S 
carries well defined symplectic structure, such that the billiard flow is the usual 
Hamiltonian flow with n commuting smooth integrals (S is the codimension two 
submanifold - the cotangent bundle of the ellipsoid, see Lazutkin [32'). Whence 
we can use the Arnold-Liouville theorem in the description of the system. Alterna- 
tively, we can consider the billiard as a discrete integrable system with the billiard 
mapping <p ■ (xfc, yt,) i-> (x^+i, yfc+i), where x^ is a sequence of the points of impact 
and Yk is the corresponding sequence of outgoing momenta |46j . However, it is not 
a simple task to explicitly describe the billiard map in the Descartes coordinates, 
as it the case when the trajectories between the impacts are the straight lines. 

The first result in this direction is performed by Fedorov 21 , who calculated 
the billiard map and found the Lax representation for the billiard system under 
the influence of the elastic force |49| . As a slight modification, in this section 
we consider the billiard system under the additional infiuence of the Rosochatius 
potential. In the derivation of the billiard mapping we use a discrete version of the 
reduction given in Section 4. 

6.2. Harmonic oscillator and ellipsoidal billiards. Consider the Jacobi 
problem (|4.3p on the 2n + 1-dimensional ellipsoid When parameter oq 

tends to zero, the Jacobi fiow transforms to the billiard problem within real 2n— 1- 
dimensional ellipsoid in C" 

(6.1) i;2"-i = {zeC"|(a-iz,z) = 1}, 

where the motion between the impacts is infiuenced by the elastic force —az (har- 
monic oscillations constrained inside the ellipsoid (|6.ip ): 

(6.2) z = p. p = — (TZ. 

Here we denoted a = diag(ai, . . . ,a„), z = (zi, . . . , z„), p = (pi, . . . ,p„). 

If (7 < 0, then all trajectories have refiections from the boundary E'^^~^, while 
for (T > 0, the initial conditions (zo,po) determining the energy h = iJ(zo,po) 
should satisfy 

/i-f- |(z,z) > e > 0, zeE^'^-K 
Consider the integral 

(6.3) J = {A-^p,p){A-^z,z) - cf{A-^z,z) 

of the Jacobi problem (|4.3p . For distinct a^, it is equal to the sum —J^i^i^^fi- 
Also, note that H/J for cr = equals to the square of the Joachimsthal integral of 
the geodesic flow [36) . 
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In the limiting process, the integral (j6.3l) multiplied by oq becomes 

(6.4) lini ao J = \{{a-^z, p) + (a^^z, p))^ + (1 - (a^^z, z)){{a-'^p, p) - cr). 

ao-i-O 4 

Let zj^^, fc e Z, be the set of impact points. By pk we denote the outgoing 
velocity at z^. From (|6.4p we get the integral 

Jfc = (a"^Zfc,PA;) + (a"^Zfc,pfc). 
of the billiard mapping within ellipsoid (j6.ip 

(6.5) (?!)(zfc,pfc) = (zfe+i,pA;+i). 

With the above notation, the mapping (|6.5I) given in |21] reads 

Zfc+i = - — [(o- - (pfe,a"^pfe))zA; + JfePfc], 

Pfe+1 = ^[(ct - (pfc,a"^Pfc))(pfe + TTfefl^^Zfe) + Jk{TTka~'^Pk - CTZfc)], 

where 

^k = \J <yJl + (o- - (Pfc, a-ipfc))^, TTfe = Jfc/(zA;+i,a"^Zfe+i). 

6.3. The Jacobi-Rosochatius billiard. The harmonic oscillations with im- 
pacts at the boundary (j6.1D have integrals 

(6.6) hj{zk,Pk) ^ -^{zk^jPk.j - Pk^jZk^j), j = l,...,n, 

and we can perform a discrete analogue of the reduction described in Section 4. 
Let us fix the values of the integrals 

hj=fij, j = l,...,n, 

where some of /ij can be equal to zero. Introduce the coordinate change 

(6.7) Zk,j ^ Xk^je'"^"- , Pk,3 ^ yk,je'^'-' + i-^e'^"- for ^,^0 

and consider the restrictions of Zkj,Pk.j to M, 

Zk,j = Xk^j e M, pfej = j/fcj e R for = 0. 
The mapping (j6.5L induces the mapping 

(6.8) $(xfc,yfc) = (xfe+i,yfc+i), Xk,y^k+i & E^'^ 
where E^^^ is the ellipsoidal component of the boundary of the domain 

£>" = {x e M" I (x, a^^x) < 1, Xj > for /^^ ^ 0, j = 1, . . . , n}. 
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Figure 1. Billiard domain for n = 2, /xi = 0, /X2 0. 



Lemma 6.1. 

Vk+lJ ■ 



where 



' (JkUkj + KkXkj) + 



-2 ' 



1 



[KkXkj + JkVkj] , 



fij ^ 0, 
fij = 0, 



t^k 
he 



Kk { —Xk,j + yk,j + 



JkC-^^"'^ 



^k 



T^k , . •"■feMj 

—yk,j - (^xk,j + 1— — - 



Xk+1,, 



Jk = 2(xfc,a Vfe), 

\Xfc Xfc 



TTfe = Jfe/(Xfe+i, a ^Xfe+i), 



6k,j = Arg 



-Kj-XkJ — JkVkj — iJk 

^k,j 



This is a billiard mapping of the billiard system within a domain £)", where 
the motion between the impacts is described by the Hamiltonian function 



1 1 " 

H{^,y) = ^(y,y) + ^(x,x) + - ^ 



j=i J 
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that is, 

(6.9) X = y, x=-(TxH 

x-^ 

where fi^/TC' = (/i?/a;f , . . . , nl/xl). 

The trajectories x(i) of (|6.9p are projections of trajectories z(t) of the harmonic 
oscillator (j6.2p (conies or degenerate conies) by the reduction (16. 7p . As above, if 
cr < 0, then all trajectories will have reflections from the boundary, while for ct > 0, 
the initial conditions (xo,po) determining the energy h = iJ(xo,yo) should satisfy 

/i+|(x,x) >e>0, xeK'-i. 

We refer to (j6.8p as the Jacobi-Rosochatius billiard mapping. After a straight- 
forward modification of Fedorov's Lax representation |21j by applying a discrete 
reduction, we obtain the following statement. 

Theorem 6.1. Suppose that the eigenvalues ai of the matrix a are distinct. Up 
to the action of the group generated by the reflections 

(6.10) {xj,yj) i-^ {sjXj,Sjyj), sj ±1, fij = 0, s^ = 1, ^ 0, 

j = 1, . . . ,n, the Jacobi-Rosochatius billiard map (16.81) is equivalent to the matrix 
equation 

(6.11) ^x. + „y. + i(A) =^x.,y.(A)£x.,y.(A)^,^\yJA) 

with 2x2 matrices depending on the parameter A 

^RioN r (\\-( Q^i^k,yk) Qxiyk,yk) + Q\i^,^) + A 
(b.iz) ^x^.y^Aj- -Qx{yk,^k) )' 



-4x„yJA) = 



/ Kk\ + JfcTTfc CrJfcA - KkTTk 

I -JfcA i^fcA 



where Q\{-,-) = (a^^-,-), oa = diag(A ~ ai, . . . , a„ - A). 

Remark 6.1. Note that the billiard can be seen as a limit of the system (|4.10p 
as ao tends to zero and /io = 0. Contrary, as the value of the integral Jk — 2(xfc, y^) 
tends to zero, the impact points Xfc approximate the trajectories of the Jacobi- 
Rosochatius problem on As in the case of the Jacobi-Rosochatius system 
(I4.10p on a symmetric ellipsoid, the Jacobi-Rosochatius billiard map (|6.8p for a 
symmetric ellipsoid (see the equation (j6.22p given below) is an example of a discrete 
system integrable in a noncommutative sense. 

6.4. The Chasles and Poncelet theorems. In what follows we shall give 
a geometric interpretation of the integrability. 

Let us recall on a well known variant of the Chasles theorem for the billiard 
system within ellipsoid 

(6.13) = {xeR"|(x,a-ix) = 1} 

without external forces, e.g., see |29L I46L I16| . Assume that the eigenvalues of 
the matrix a are different. Let x^ be a generic sequence of impact points. Then 
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the sequence of lines x^x^+i is simultaneously tangent to the same set of quadrics 
Qi)i, ■ • • , Qj7„_i from the confocal family 

n 2 

(6.14) E^T^-l^ 

2—1 

Namely, the set of lines p = p(x, y) = {x + sy | s e R} that are tangent to the 
quadric from the confocal family (|6.14p are given by the equation [36) 

(6.15) $x,y(?7) - (gr,(x, x) + l)g^(y, y) - (x, y) = 0. 

On the other hand, from (I6.1ip where we set cr = and = 0, $xfc,yfc(A) = 
det£xfc,yfc(A) is an integral of the system. Therefore, if 77 is a zero of <&Xfc,yfc(A), 
then the lines pk — p(xfc,yfe) = x^x^+i, fc S Z are simultaneously tangent to Q,,. 
For a generic trajectory (xfc,yfc), we have n — 1 distinct solutions 771, . . . ,77,1-1 of 
the equation <l?xfc,yfc(A) = 0. 

Moreover, suppose that the trajectory x^ is periodic. Then any billiard trajec- 
tory which shares the same caustic quadrics is also periodic, with the same period 
{the Poncelet theorem in R" [59ll9l [T5l[T6] ) . 

An analytical condition on caustics Qr^i, ■ . ■ , Q?7„_i for periodic billiard trajec- 
tories is derived by Dragovic and Radnovic, generalizing classical Cayley's condition 
for n = 2 |15L 116) . The geometry of the lines common to the confocal quadrics 
is further studied in |30L 116) . while Chasles's-type theorems for several natural 
mechanical systems are given in |33L I19L 122] . 

Let d be the number of indexes i for which fii ^ 0. 

Theorem 6.2. Assume that the eigenvalues of the matrix a are different. Let 
(xi;,yi;) be a generic trajectory of the Jacobi-Rosochatius billiard map (16. 8p . 

(i) Assume a ^ (respectively, a = 0). The trajectories 

(6.16) h: xfc(t), teR 

of the Jacobi-Rosochatius system (|6.9p with initial conditions (xfc(O), yfc(O)) = 
{x.k,yk) are simultaneously tangent to the quadrics 

(6.17) Qrii , ■ ■ ■ , Qrin+d {respectively^ 7 • • ■ i Qrin-i+d) 

from the confocal family (I6.14p . where rji, . . . ,r/n+d (respectively, r/i, . . . ,r]n-i+d) 
are solutions of the equation det Cx^^y^{X) = 0. 

(ii) Suppose that the trajectory (xfc,y/j) is periodic. Then any billiard trajectory 
which .shares the same caustic quadrics (j6.17p is also periodic with the same period. 

Proof, (i) Firstly, we take a = 0. For a sake of simplicity, assume /ii, . . . , /i^ 7^ 
0, ^id+l = ■ ■ ■ = fin = 0- 

Consider the billiard system without external forces within a symmetric ellip- 
soid 1 in 

R"+^ - C^(zi, . ..,zd)x R"-'^(xd+i, . . . , Xd), 
where we take a = (oi, ai, . . . , Od, Od, Od+i, 0^+2, . . ■ , ««), Oi ^ aj, i ^ j. Let 

(6.18) (Xfe,yfc) = (2l,fc, . . . , Zd^k, Xd+l,k: ■ ■ ■ ,Xn,k,Pl,k, ■ ■ ■ ,Pd,k,yd+l, ■ ■ ■ ,yn,k) 
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be a generic billiard trajectory with the values of the integrals 

(6.19) hj -^{zk^jPk.j - Pk.jZk^j) = tij, j = l,...,d. 

From Lemma (|6.2I) below, the lines pk = XfcXfe+i C x R"~'^ determined by 
the billiard trajectory (|6.18p are tangent to the N = n + d — 1 quadrics 

(6.20) Q,. : + ■ • • + + +■■■ + - 1, 

where 77; are zeros of the corresponding polynomial (j6.23l) . 

By applying the reduction (16. 7p in the variables zi , . . . , we get the billiard 
trajectory (xfc, y^) of the Rosochatius billiard map (|6.8I) with a — 0, fii, . . . , 7^ 
0, Hd+i = • • • = = 0. At the same time, det -Cx^.y^ {r]i) = 0, i ^ 1, . . . , n + d — I 
and the lines pk project to the curves (|6.16p tangent to the quadrics (16. 17^ . 

In the other direction, let (xfc,yfc) be a trajectory of the Rosochatius billiard 
map (16. 8p with a = 0, /Lti, . . . 7^ 0, /id+i = • • • = /i„ = 0. Then we can lift 
(xfejYfe) to the 5'0(2)'^-invariant family of billiard trajectories (|6.18p of the billiard 
system without external forces within £'"+''-1 satisfying (j6.19p . 

The case cr ^ 0,^i = • • • = pn = is proved by Fedorov |21l . For a ^ 0, 
instead of the equation (j6.15p . which characterizes tangent lines to the quadric Q,,, 
we use the equation 

(6.21) y (,7) - (Q,(x, x) + l)(Q,(y, y) + a) - Q^x, y) = 0. 
A conic / = {x{t) \ t eM.} associated to a solution of the equation 

Mt) ^ y, y{t) -crx 

with the initial condition (xo,yo) is tangent to Q,, if and only if "I>Jy(77) — 
(Proposition 1 in [21 ). The statement can be applied for matrices a with multiple 
eigenvalues and repeating the arguments given above, item (i) follows. 

(ii) Suppose that the trajectory (x^, y^) with the energy H — his periodic. The 
Jacobi-Rosochatius billiard (|6.8p can be considered as a usual integrable Hamilton- 
ian system on T*D/r \ E (we take only trajectories which are not tangent to the 
boundary E^~^, that is Jk — 2{xk,yk) 7^ 0), see |32| . The associated continuous 
billiard trajectory (x(t),y(i)) is periodic as well. 

The parameters of caustics 771, ... , rjn+d (respectively, 771, ... , rjn+d-i) corre- 
spond to the unique values ci, . . . ,c„, ci + • • • + c„ = 2h oi the integrals (j4.12p 
(where we set i,j — 1, . . . ,n) determined from the equations 

n n 2 

The curves (j6.16l) tangent to the caustics (|6.17l) are associated to billiard tra- 
jectories which belong to the Lagrangian tori that are components of the invariant 
level set 

Me,,...,c„ C T*D/r \ S : = ci, ...,/„= c„. 
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According to the Arnold-Liuville theorem, if a trajectory on a regular La- 
grangian torus is periodic, all trajectories on the same torus are also periodic with 
the same period. Since all components of Mcj_..._c„ are related by the reflections 
(|6.10p . it is clear that if one Lagrangian torus is periodic, so are the others one. □ 

We note that there is a natural generalization of the presented results to the 
ellipsoidal billiards on spheres and hyperbolic spaces |29L I47L I14L I42j . 

Lemma 6.2. [The Chasles theorem for a symmetric ellipsoid] Consider a billiard 
within a symmetric ellipsoid E"~^ defined by 

(6.22) a, i G /i {1, . . . , fej, . . . , 

= ar, i e /,■ = {fci + • • • + fcr-1 + 1, . . . , fco + • • • + fcr-1 + , 

Ui ^ aj,i ^ j , ki + ■ ■ ■ + kr — n. Let be a generic sequence of impact points. 
Then the sequence of lines XfeX^+i is simultaneously tangent to the same set of 
quadrics Q^-^ , ■ ■ ■ , Qijjv from the confocal family (|6.14p . where N = Si + ■ ■ ■ + Sr — 1, 
6s — 2 for ks — \Is\ > 2, and Ss = 1 for k^ = \Is\ = 10 Further, due to the 
SO{ki) X • • • X SO{kr) symmetry of the system, if Xk is a billiard trajectory, so is 
i?(xfc), R e SO{ki) X ••• X SO{kr), and the lines R{'x.kXk+i) are tangent to the 
same set of N quadrics. 

Proof. Note that the description (|6.15p of the tangent lines to the quadric 
does hold for the matrices a with multiple eigenvalues (|6.22p . According to (|5.2I) . 
the number of quadrics tangent to a generic line p — p(x, y) equals to the number 
of zeros of the polynomial 

*x.y(A) - (A-ai)^i •••(A-a,)^'-det£x,y(A) 

(6.23) = E ( - "^)'^"' - "^)'"/^ + 

where fs, Ps are given by (|5.3p . (|5.4I) (where we set cr = 0, /i^ = 0). Therefore, a 
generic line p = p(x, y) is tangent to N — 6i + ■■■ + Sr — 1 quadrics from the family 

(EH- 

Consider a generic billiard trajectory (xfe,yfc) within a symmetric ellipsoid 
(I6.13|) . (|6.22p . As above, if the line pk = p{xk,yk) — x^x^+i is tangent to the 
quadric defined by the equation (|6.14l) . i.e, ^Xfc,yfc('?) = 0, then the lines pt, 
k gZ are also tangent to Q^. □ 




Note that N = r + p — 1, where r + p is the dimension of invariant isotropic tori of the 
corresponding geodesic flow on the symmetric elUpsoid HI. Ill in R"+^ with semi-axis ^/ao ~ and 
^/oi, . . . , Y^o^ given by 116.221 1 (Theorem [5T2j . 



24 



JOVANOVIC 



7. Hierarchy of the Lelx representations 

By taking the limit ao — and assuming /ip — 0, from Theorem 14.11 we can 
write down the Lax representation 

£x,y(A) = [£x,y(A),^x,y(A)], 

QA(x,y) Qx{y,y) + Qx{^.^) + o 



for the Jacobi-Rosochatius system on E" (|6.9p . i.e, for the trajectories of the billiard 
system between the impacts. Here we used 



(x,a-2x) + ^(l-(x,a-ix)) 

as ao —J' 0, where x — (a;o,x) and y — {yo,y) satisfy the constraints (jl.4l) . 

This Lax representation, for fxi — ■■■ — fin, is equivalent to the Lax repre- 
sentation for the harmonic oscillator given in [17] . where, by induction, the Lax 
representations for natural mechanical systems with polynomial potentials separa- 
ble in elliptic coordinates in M" are given (see also [2]). In the same way, by using 
Theorem 14. 11 we can give the 2 x 2-Lax representations for the separable potential 
perturbations of the Jacobi-Rosochatius system on £'". 

7.1. Separable potentials in Recall, a potential V{x) is separable in 

the elliptic coordinates Aq < ap < Ai < ai < • • • < A„ < a„ defined by (|1.7p if and 
only if it is a solution of the Bertrand-Darboux equations 

( d d \ ( "^^ dV \ 



(Benenti |4], see also Marshall and Wojciechowski |34| ). The solutions of (|7.ip can 

dv 



be written in the form V{x) ~ Xi-S^, where V(a;) are solutions of 



(see |51| ). Then, a complete set of commuting integrals is given by 
(7.2) /,^y2^^I^!^l^^ + 2F,(a:), z = 0,l,...,n, 

where F,{x) ^ x,§§- 

Polynomial potentials are described in [501 151L I27j . Basic homogeneous poly- 
nomial solutions V'^^^ of degree 2k of the equations (j7.ip in elliptic coordinates 
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reads 

....-f^t^n^ 

and V^'^'' , Fq''^ , • ■ • , Fn''^ , fc € N satisfy the system of the recurrence relations 
(7.3) Ft+'^ = a,F^'^ ~ x^V^'l F^'^=xf, V^''\x) ^ F^"^ + ■ ■ ■ + Fi''^ 
(we use the notation given by Zaitsev I51j ). For example, 

V^^'> = {x, x) (the Hook potential), ^ = x"^, 

= {Ax, x) - V'^^Xx; x) (the Garnier potential), F^;^^ = x1{a, - V^^'), 

The rational and the Laurent polynomial solutions of (j7.ip are given in [50L 
I27L I12L I13j . The basis for degrees —2 and —4 are given by 



= ^ (the Rosochatius potentials). 



s = 0, . . . , n. 



Qi - asX% • flj - fls a;| \ ^ - a., 

7.2. Natural mechanical systems on ellipsoids with separable poten- 
tials. Consider the motion of a material point on an ellipsoid (jl.ll) under the 
influence of the potential 

(7.4) v{x) = V+{x) + -Y,^, F+(x) = -^a,yW(x), 

1=0 * k=l 

where ak are real parameters. The equations of motion are 

(y.. . ■ {A-^y,y)-{WV{x),A-^x) 

7.5 x^y, y = -— ^ ^ A x-VVix). 

[A ■'x,x) 

As a straightforward generalization of Theorem 14.11 by using the constraints 
(|1.4p . the recurrence relations (17. 3p . and the identities 

fe-i 

yl;,^''' = (A - A)-^A'' = A'-'A"^ - J2 A''"'^', e N 

we get. 
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Theorem 7.1. Suppose that the eigenvalues ai of the matrix A are distinct. 



Up to the action of the group generated by the reflections 

[xi,yi)^ — > {siXi.Siyi), = ±1, i = 0, 
the system (|7.5p is equivalent to the matrix equation 
(7.6) t{\) = [C{\),A{X)] 

with 2x2 matrices depending on the parameter A 

C(x)=( <i>^i^^y) '?A(2/,y) + gA(f,f) + A(a;,A)\ 
V^l - '?A(a;,x) -qxiy.x) J' 

_ 1 f i((VF+(x), A-ix) - {A-^y,y) (^-^f , f )) ~ A)^ 
' {A-^x,x) {{A-^x,x) , 

w/iere 

A(a;, A) = o-iAi(a::, A) H h (t„A™(x, A), 

Afe(x, A) = A*^-i - A'=-V(i) - A'^-^yia) aV^'^-^) _ p.(fe~i)^ 

ri(x, A) = (A^^a;, a;) (f7ifii(a::, A) H h CTmr2m(a;, A)) , 

and Qkix,X) are determined from the equations 



2nk{l + qxix,x))^2Akix,X) + {A-'x,VV^''\x)), k = l,.. 
For example, 

Ai = l, = 
A2 = A - (a;, x), n2^X-2{x,x), 
A3 = A^ - X{x,x) - {Ax,x) + {x,x)'^, 
ris^ X^ ~ 2X{x,x) - 2{Ax,x) +3{x,xf. 

Now, the integrals (|7.2p and the Lax representation are related by 



, m. 



det(£(A)) = il + qxix,x))(^qx{y,y) + qx(^^,'-^) +A{x,X)'^ -qxix,y) 

2 



//fc I L p lb 



A — fli ' ^ (A — ai)2 

fe=l 1=0 i=0 ^ ^ 

By taking ~ a I2Y^'^^ ~ a/2{Ax,x) — a/2{x,x)^ and assuming /ig = 0, in 
the limit oq — ^ , we get the Lax representation for a natural mechanical system 
in R" under the influence of the the Garnier potential V = (T/2(ax, x) — cr/2(x, x)^ 
obtained by Antonowicz and Rauch-Wojciechowski [2] (see also |41| ). Note that, 
similarly as in Eilbecktt, Enol'skii, Kuznctsov and Tsiganov [17] . one can consider 
the problem within a framework of r-matrix method. 

Finally note that the polynomials V^''\ as well as the Lax representation, 
are well defined for a symmetric ellipsoid (|l.ip . (15. ip . Therefore, repeating the 
construction presented in Section 5, we obtain. 
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Corollary 7.1. The system (17. 5p on a symmetric ellipsoid p.ip . (|5.ip is 
completely integrahle in a non- commutative sense by means of integrals (j5.5p anrf 



Generic trajectories take place over r + p-dimensional invariant isotropic tori, 
spanned by the Hamiltonian vector fields ,Xp^. Also, the problem is Liouville 

integrable by means of integrals fa and (|5.15p . 
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